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Abstract. We study the geometry associated to the Grusin operator 

G = A X + \x\ 2 dl on K™ x M„, 

to obtain heat kernel estimates for this operator. The main work is to find 
the shortest geodesies connecting two given points in R n+1 . This gives the 
Carnot-Caratheodory distance dec, associated to this operator. The main 
result in the second part is to give Gaussian bounds for the heat kernel Kt 
in terms of the Carnot-Caratheodory distance. In particular we obtain the 
following estimate 

MC,r,)\ < Ct-9" 1 min ( 1 + 1 + Y e - ^cciCvf 

\ \x + £\ At J 

for all C = (a, Ml), V = (£>«) 6 where a = max{f - 1,0}. Here the 

homogeneous dimension is q = n+2, so that \ — ^- = ^-?r- This shows that our 
result for n > 2 corresponds with the result on the Heisenberg group, which 
was given by Beals, Gaveau, Greiner in [T|. 
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1. Introduction 

The purpose of this article is to study the geometric properties of the Grusin 
operator 

G = A x + \x\ 2 d% onW>xR u 
to give estimates for the heat kernel of this operator. One may write G as 

71 

with the smooth vector fields 

Xj := d Xj , Uj := xjd u , j = 1, . . . , n. 

Note that G is hypoelliptic, since G satisfies the Hormander condition. We are 
interested in the Carnot-Caratheodory distance associated to these vector fields, 
i.e. the length of a minimizing horizontal curve. A horizontal curve is a curve, 
whose tangents are a linear combination of the vector fields Xj, Uj, j — l,...,n. 
Depending on the start and end point of the curve, there are a different number of 
locally minimizing curves, i.e. geodesies. In addition to the results of Calin, Chang, 
Greiner and Kannai in [3], we can compare their lengths to give explicit formulas 
for the Carnot-Caratheodory distance. The geodesies starting in £ R n+1 are very 
similar to those on the Heisenberg group (see [Tj). But comparing the lengths of 
geodesies starting in a generic point (x, 0) eR"xK gets more complicated. Another 
difference occurs when studying curves connecting (x, 0) and (±x, it), where on gets 
two different types of geodesies (see Theorem (|2.4p and l|2.5p ). 
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These results on the Carnot-Caratheodory distance allow us to give some point- 
wise estimates of the heat kernels. 

Due to the fact that the partial Fourier transform of G in u leads to the rescaled 
Hermite operators, it is easy to calculate the heat kernel using Mehler's formula, 
up to the partial Fourier transform. One may also observe that G is translation 
invariant in the u-variable. This allows us to write the heat kernel in the form 
Kt(x, £,it), x,£ £ R™, uel. So the solution to the heat equation 

Ud t ~G)u =0 
\«(0,-) = / 

is given by 

u(t,x)= [ f K t (x,£,u-\)f(£,\)d\dS. 



We are able to give global estimates for the kernel of the form 

\K t (x,£,u)\ < F(x,^u,t)e-^ dcc({xfi) ^- u))2 , 
where F is a function with polynomial growth. It turns out that 



2 \ a 



F(x,£,u,t) < Ct * mm I 1 + -r^-j-^-, H — I 

with a = max{| — 1,0} and C > is a constant independent of x, £,u,t. In 
particular our main result is the following 

Theorem. For £ = (x, 0), r/ = (£,u) € R Tl+1 we have 

(1.1) \K t (x,Z,u)\ < *-f-min (l + ^f^,l + d -2£^J e -^^\ 

with a = max — 1,0). 

For a better understanding of the exponent a, consider the dilation 

6 r (x,u) := (rx,r 2 u), x G R", u e R, 

for r > 0. Then G is homogeneous with degree 2, and also the Carnot-Caratheodory 
distance is homogeneous with degree 2, i.e. 



JTI + I 



In this setting the homogeneous dimension is q = n + 2, and § — 1 = f + 1 = §■ 
This shows now that our result exactly coincides, in the case n > 2, with those on 
the Heisenberg group obtained by Beals, Gaveau, Greiner in [l]. 

There are some Gaussian-type estimates for heat kernels with greater generality. 
Sikora, for example, shows in [8] 

\Kt(x,€,u)\<t =11+ — I e 

The methods that are used are completely different. The explicit formulas in our 
situation allow us to give a smaller exponent of the polynomial factor, i.e. 
instead of . 
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2. The sub-Riemannian geometry associated to the Grusin operator 

Let M be a connected C°° -Manifold with smooth real vector fields Xx, . ■ ■ , X m . 
For x G M and v <E T X M define: 

!m m 
3=1 3=1 

where we use the convention, that inf := oo, i.e. if v is not in the linear span 
of Xi, . . . ,X m , then we set \\v\\ x := oo. A horizontal curve 7 : [0,1] — > M is 
an absolutely continuous curve such that 7(f) € span {.Xi (7(f)), . . . , X m {^{t))} for 
almost every t € [0, 1]. We define the length of a horizontal curve by 

1 

(2.2) m--= J \m%( t) dt. 



Then the Carnot-Caratheodory distance of two points p, q € M associated to 
X\ , . . . JT m is defined by 

(2.3) d cc (p, 9) := inf L( 7 ), 

where the infimum is taken over all horizontal curves 7 connecting p and g, i.e. 

7(0) =p, 7(1) = 

It is well known by Chow's theorem (see [1]) that, if X\ t . . . , X m and their brack- 
ets span the tangent space T X M at every point x of M, then any two points can be 
joined by such a curve, so that dcc(p, q) < 00, for any p,q € M. 

Now let M = and 

= 9 Kj , Uj = Xjd u , j = 1, • • • , n. 
Then the Grusin operator G = A x + \x\ 2 reads in terms of these vector fields 

n 

The vector fields X±, . . . , X n , U\, . . . , U n span the tangent space everywhere except 
along the line \x\ — 0. But, since [Xj,Uj] = d u , j — l,...,n, all conditions are 
fulfilled to connect any two points by an absolutely continuous curve with finite 
length. If 7 : [0, 1] — > R™ +1 is an absolutely continuous curve, one has 

n 

7 = in+id u + 

3=1 
n 

= jn+idu + Y ii X i 

3=1 

almost everywhere. Since the vector fields Uj, j = 1, . . . , n are not linear indepen- 
dent, we may write 7„+i9 u as a linear combination of U±, . . . , U n in different ways. 
So 

n In 

Il7(*)ll?(t) = E ^i(*) 2 + inf E a ) ' 7«+i(*) = E a i7i(*) 
3=1 [3=1 3=1 

To minimize this convex functional, one can use the method of Langrange multiplier, 
which gives 

I ■ /-All 2 7n+l W , • M2 
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Now let 7 = (7(1), 7(2)), with 7 (1) = (71, . . . ,7„), 7 (2) = 7«+i, then 

I, • /,x|,2 7(2) it) 2 I . , ,|2 

h(*)ll 7 (i) = 1 — [2 + |7(i)(*)| • 

|7(i)(*)| 

To calculate the Carnot-Caratheodory distance, one can minimize the "energy" 
integral 

1 1 / • 2 \ 

(2.4) J ii7(*)ii; ( t)*= / (f^^+i^i 2 )^ 

^ |7(i)l /| / 

This leads to the Euler-Lagrange equations, which gives us local minimizing curves, 
i.e. geodesies: 

(2.5a) 7(1) + 7^37(1) = 

|7(i)| 

1 I 

Solving these equations yields the following result (see [SI El [7]). 

All geodesies 7 starting in (xi,wi) el"xl are given by r y b ' c — (7(1)) 7(2) J ; 

c 

- sinfbt) + X\ coa(bt) 
b 

^/t_sin(2W)\ ^ 
b \2 46 J b 
a, ft sm(2bt)\ 

with parameters 5g R, 6^0, c <E R™, and 

(2.7) 7°' C W = (ct + xi.m) 

which is the limiting case 6^0. 
The length is given by 



[o, 1] -> 


R" x K, where 


(2.6a) 


7?i)W 


(2.6b) 


7( b 2 )(*) 



(2.8) L(7 b ' c ) = ^/icf + NV. 

To calculate the Carnot-Caratheodory distance, we are interested in the shortest 
geodesic joining two given points (xi, Ui), (a;, u) € K" x R. 

In our further analysis we will assume that u± = 0, since 7 is a geodesic con- 
necting (xi,0), (x,u), if and only if 7 + (0, u\) is a geodesic connecting {x\,u\), 
(xi, u + u\); with the same length. 

We will also assume that u > 0, since 7,^ = 7^)'° an d 7( & 2 ) — — 7(2)' C - 

Furthermore observe, that given 6 6 1\ 7rZ* the parameter c € R™ is uniquely 
determined by c = rj^r (x — x\ cos 6), due to the boundary conditions. We will 
denote this geodesic by 7 b for short. 

Given b € 7rZ* (note that this means x = ±xi), the parameter c € R™ is 

determined by the whole sphere |c| = \J2bu — \xi\ 2 b 2 . But any c 6 R™ with 



|c| = \J2bu — |xi| 2 6 2 will give the same length. 

Theorem 2.1. Given (xi,0),(x,0) € R™ x R, the only geodesic connecting them 
is given by 

(2.9) 7°(i) = (x 1 +*(x-x 1 ),0) 
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and its length, which is equal to the Carnot-Caratheodory distance, is 
(2.10) L( 1 °) = d C c((xi,0),(x,0)) = \x- Xl \. 

Proof. We have to prove that b = is the only parameter that allows connecting 
(xt,0) and (x,0). But from the Euler-Lagrange equation, we know that 7( 2 ) = 

I 1 2 

b 7(i) ■ That means 7( 2 ) (t) > or 7(2) (i) < for all t > 0. But since we want 
7(2) (0) = and 7(2) (1) = 0, we conclude that 7(2) = 0, so that 6 = 0, hence 

C = X — X\. □ 

Theorem 2.2. Given u > 0, there are infinitely many geodesies -y m7r < c >" ; meN, 
connecting (0,0) € E n x M and (0,u) G R" x M jjscn by flU} toft/i c m € R", 
|c m | = \J2mwu. Their lengths are 



(2.11) 



L( 7 m7r > c 



V2 



m/KU. 



In particular the shortest geodesies in this case, which give the Carnot-Caratheodory 
distance, are ~f^ Cl , \c\\ = \j2nu; so 

(2.12) dac((0,0),(0,u))=V27rtt. 

Proof. The boundary conditions 7^(1) = and 7p)(l) — u are equivalent to 

sm(2bY 



c c 
= - sin 6 and u = — — 
b 2b 



1 - 



2b 



Since u ^ 0, we have c 7^ 0. But then there exists m £ N with 6 = m7r, so that 



2mnu. This gives the claim. 



□ 






IOtt 




Figure 1. Geodesies starting at (0,0) and ending in (0, u) with 
different parameters. 



Now we consider the case where we want to connect (a;i,0) and (x,u) with 
x 7^ ±xi, u > 0. The first lemma is just a rewritten formulation of the boundary 
conditions 7 fc ' c (0) = (afi,0) and 7 b,c (l) = (x,u). 

Lemma 2.1. Suppose (aij.,0), (x,u) £ R n x R, x 7^ ±£1, it > 0. TTien 7 b ' c given 
fry l|2.6|) is a geodesic connecting these two points, iff 

b 

rx — Xibcotb 



sin b 



and b is any solution of 
2u 

(2.13) 



! 2 + |d 2 sin 2 b 



cot b 



2x ■ x\ 1 — b cot b 
|xi| 2 + |x| 2 sinb 



HEAT KERNEL ESTIMATES FOR THE GRUSIN OPERATOR 



6 



The square of the length of r y b = r y b,c is then given by 

b 2 / 

(2.14a) i 2 (7 fc - c ) = — o- [\xi\ 2 + \x\ 2 - 2x -xicosb 

sin b V 

(2.14b) = 2bu+ [\xi\ 2 + \x\ 2 ^bcotb-2x-x 1 -J^. 

Proof. To find those geodesies connecting two points (#1,0) and (x,u), we have to 
find b ^ and c € R n , such that 7({)(1) = x and 7( 2 )(1) — u - Since b G 7rZ means 
a; = ±xi, we have that b $ 7rZ. Now 7^(1) == u gives: 

|c| 2 / sin(2b)\ „%i-c . 9/ ,. . ,9 , / sin(26) 
2. = if ^1 - -L2 J + 2— sm» + N b(l + ^> 

(2.15) =\ (|c| 2 + |x x | 2 - ^ (|c| 2 - l^l 2 ^ 2 ) + 2^ sin 2 6 

6 



And 7,^(1) = x gives 



(2.16) -x — c + x\bcotb, 

smb 



so that 



b 2 



— l_|a;| 2 = |c| z + 2c-xi6cot6+ \xi\ z b 2 cot' 2 b 
sin 6 

(2.17) = \c\ 2 + 2c-x 1 bcotb-\x 1 \ 2 b 2 + |a:i| 2 62 

sin & 

and 



I i2 & 2 , ,2 6 2 COt6 2 2 2 

C = n— \x\ — 2X\ ■ X ; — ; h \Xl\ b cot 

sin^ b sin 6 

b 2 /, ,9 , ,5\ - & 2 cot6 

sin 2 b 
b 2 

(2.18) = — 5- (|x| 2 + |xi| 2 -2x-xiCOS&) - |xi| 2 & 2 

sin 6 V / 



(\xf + \x 1 \ 2 )-2x 1 -x b -^-\x 1 \ 2 b 2 
\ I sin 



sin b 

We can continue with (|2.15p by using equation l|2.17p 

1 A |2 , 1 ,2, 2 \ A |2 , ,2\ , „C-Xi 2 , , n c ' x l 



2u = i (|c| 2 + |x!| 2 6 2 ) - cot b (\x\ 2 - l^il 2 ) + 2°—^- cos 2 6 + 2^ sin 2 b 
= i (|c| 2 + |x x | 2 b 2 ) - cotb (|x| 2 - \ Xl \ 2 ) + 2- 



C • Xl 



1 /, ,2 . I ,2, 2 \ , , /, ,2 ,2\ C- Xi + |xi| 2 6cot6 

= - (Jc| +\xi\ b 2 j~cotb[\x\ +|xi|J+2 . 

Now use equation (f2TT6]) and (f27T8j) 

x • Xi 



2u = 



±(\c\ 2 + \x x \ 2 b 2 ) -cot6(|x| 2 + |x 1 | 2 ) +2- gii] , 

b A /, ,j . ,2\ 1 — bcotb 

-5 coto X + xi + 2x • xi : — . 

, sin b J V / sin 

And the square of the length is given by 

L 2 {^ c ) = -^5- f|x| 2 + |x 1 | 2 -2x-xiCOsb 
sin 6 V 

= 26u + bcotb f |x| 2 + Ixil 2 ) — 2 x • x\. 

\ J sin 6 



□ 
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Definition. For b e R\ ttZ*, -1 < a < 1 define: 

n b 1 — 6 cot b 

(2.19) mO) : = MM) := — cot6 + a- 



sin 2 6 sin 6 

sin 2 6 

In the following we will often use the abbreviation 



(2.20) 1(b) := l(b,a) := — — (1 -acosb). 



R 



= V 1^1 +M and a 

Fi| + F 



Next we will study the functions /x and /2, /t, which will be introduced in the 
next lemma. (Note that fx in this section does not mean Fourier transform!) 

Lemma 2.2. For b € K \ 7rZ*, < a < 1, /i is a convex combination, i.e. 

(2.21) Li(b) = (1 - a)/x + a/t ' " 
and for — 1 < a < 0: 

(2.22) ^0) = (1 + a)M( fo ) - a M 
wftere 

(2.23) /2(6) := — \ cot 6, be K\ttZ*, 

sin 6 

(2.24) /t(6) := — b — +tan&, 6 e R\ I ( k+ I ) tt : fc € Z 

cos^ 6 ^ V 2/ 

Proo/. For 6 e R \ ttZ* : 

1 — b cot 6 H — r^-r — cos 6 
M (6) = (1 - a)H(b) + a 

sin 2 ^ Cl 4- -M 
= (1 - a)A(6) + 2a 2^+sin 6 J 



and 



sin b 

(l-a)£(&) + ^— I^ + tan^ 
(1 - a)jl{b) + aLi [\ 



s , v cot 6 — 1 + — cos 6 
/i(6) = (1 + a)£(&) - a 



sin 6 



cos 



2 b ( b 



2 (sinfc "0 



= (l + a)/2(&) - 2a 



sin 6 

- (1 + a)£(&) - a [ - cot ^ 



(1 + a)fi(b) — afi 



□ 



Remark. Observe that all these functions fi, jl, ft are odd. The following lemmata 
will show that /i|[o,oo) > 0j AI[o.oc) > an d therefore At|[o,oo) > 0- Since we assumed 
that u > 0, we only need information of these functions for b > 0. 
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7T 2-7T 37T 47T TT 2lT 3tT All 

(a) Plotting of jX(b) (gray), p.(b/2) (gray) (b) Plotting of p.(b) (gray), ju(6/2) (gray) 

and /i(b) (black) with parameter a = and fi(b) (black) with parameter a = 

0.7. -0.7. 



The function p,, which also appears in the study of geodesies on the Heisenberg 
group, where intensively studied by Beals, Gaveau, Greiner in [l]. We will use the 
following result: 

Lemma 2.3. The function p is a monotone increasing diffeomorphism of the in- 
terval [0, 7r) onto [0, oo). On each interval (mir, (m + l)ir) , m € N, p has a unique 
critical point b m , which is a minimum. b m is implicitly given as the solution of the 
equation 1 — bcotb = in the corresponding interval. On this interval, p decreases 
strictly from +oo to p(b m ) and then increases strictly to +oo. Moreover, 

(2.25) p(b m ) +7T < ft(b m+1 ), 
and 

(2.26) p,(b m ) > rmr. 

The function ft has been studied by Calin, Chang, Greiner, Kannai in [3]: 

Lemma 2.4. The function ft is a monotone increasing diffeomorphism of [0, 5) 
onto [0, oo). On each interval [mm + \, (m + 1)tt + f) , m € No, fx has an unique 
critical point b m , which is a minimum. b m is implicitly given as the solution of the 
equation 1 + 6 tan b = in the corresponding interval. On this interval, pi decreases 
strictly from +oo to p(b m ) and then increases strictly to +oo. Moreover, 

(2.27) A(6m) > T (m + i V 

Remark. Moreover p,"(b) > 0, 6 £ 7rN and /i"(6) > 0, 6 - | ^ N , i.e. /2 is 
strictly convex on each interval (mir, (m + l)7r) , m 6 N and /I is strictly convex 
on each interval (mir + |, (m + l)7r + |) , m 6 No- This is another result of [l, 
3j. Hence \x, as a convex combination, is also strictly convex on each interval 
(mir, (m + l)w) , m € N. 

Combining the previous two lemmata we get the following result for the function 
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Lemma 2.5. Let — 1 < a < 1. The function [i is monotone increasing on [0,7r) 
onto [0,oo). On each interval (rrm, (m + 1)tt), m G N, /x has an unique critical 
point b rn , which is a minimum. For a > and m € No 

(2.28) (2m + 1)tt < 6 2 m+i < hm+i < 2b m < (2m + 2)tt, 
(2m + 2)7r < &2m+2 < &2m+2 < (2m + 3)7r; 

anrf /or o < and m G N 

(2.29) 2m7r < b 2m < b 2m < 2b m < (2m + 1)tt, 
(2m + l)7r < 6 2m+ i < o 2m+ i < (2m + 2)ir. 

On each interval (mir, (m + l)7r), m € N, /i decreases strictly from +oo to [i(b m ) 
and then increases strictly to +oo. Moreover for all m € N, 

(2.30) /i(&m) > 

(2.31) l{b m ) - b m n{b m ) = l-aS(b m ), 
where 

(2.32) <5(6) := cos 6+ ^sin6. 

Proof. It is easy to see that mir < b m < mir + § < b m < (m + 1)tt, see figure [2| 




Figure 2. The graphs of the functions tan (black) and 
(gray) and the line x>—>x. 



cot 



Since 



1 — 6 cot 6 



> 




6 6 
- - cot - 
2 2 

6 6 ((2m + 1)tt, (2m + 2)tt) 
b 6 e (2m7r, (2m + 1)tt) 
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and x i— > 1 — x cot x (gray) . 



and 6 2 m+i, 26 m 6 ((2m + l)7r, (2m + 2)7r) it follows that &2m+i < 26 m ; also, since 
62m, 26 m G (2m7r, (2m + l)7r), it follows that 62m < 26 m . 

If a > 0, fi' is a convex combination of p! and O n the interval 

((2m + l)7r, (2m + 2)71") both functions are monotone increasing with zero b2 m +i 
of p,' lying left of 2b m , so we have 

M' (&2m+l) = I A' ( < and M' (2&m) = (1 - 0)A' (26m) > 0. 

Hence p! has a zero b m in between: b2 m +i < b m < 2b m . Due to the fact that fi is 
strictly convex, this critical point is unique in the interval ((2m + 1)tt, (2m + 2)ir), 
and also a minimum. 

In the interval ((2m + 2)7r, (2m + 3)tt) still both functions are monotone increas- 
ing, but only \j! has a zero at &2m+2- |/*' (|) ^ s positive on this interval and since 
p!(b) tends to —00 as 5 fc>(2m+2)^ ^ ^ there is a zero &2m+2 of // in this 

interval ((2m + 2)tt, (2m + 3)7r) with (2m + 2)7r < &2m+2 < &2m+2- This critical 
point is unique and a minimum, since fi is convex. 

The case a < is similar, if we use \ji (|) instead of i/t (|) and the fact that 
in the interval (2m7r, (2m + l)7r): 62m < 26 m are zeros of these two functions and 
\p (I) has no zero in ((2m + l)7r, (2m + 2)77). 

To obtain the lower bound (|2.3Q[) on fi, we consider the following cases: 

1. case a > and m even: 

Since *f e ((2^2 + ±) tt, (f + |) tt), we have 

M(6 m ) = (1 - o)/2(6 m ) + a/t 
m - 

> (1 — ajmn + a — - — 7r 
m — 1 
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2. case a > and to odd: 

Since h -f e ((fflfi + \) n, ^tt), we have 



to 

l^{Omj > (1 — ajmn + a-^ft 

to to—I 
^ 2"^ > 



3. case a < and to even 

2 c I 2 ' 2 



Since % G (f tt, ^tt), we have 



M(&m) = (1 + a)(l(b m ) - a/2 ( 



m 

> (1 + ajmn — a "^" 7r 
m m — 1 

> — 7T > 7T. 

- 2 2 



4. case a < and to odd 

b m r ( m-1 raj 
2 c I 2 " ' 2 



Since V G l 2 ^ 71 ", ^T 1 ^ we have 



m — 1 

A*(6 m ) > (1 + ajrriTT — a — - — n 
to — 1 

At last, to establish the equation (|2.3ip . we need some more information of b m , i.e. 
the critical points of fi. First observe, that for any 6 6 K \ 7rZ*: 

(2.33) 6/j(6) + 6cot6 b —a = l{b) 

sm 6 

and 

d /, , 6 \ , b 1 - 6cot6 
o cot b -a = cot 6 ^ a- 



So 



d6 V sin b J sin 2 6 sin b 



l'{b) = fi(b)+bfx'(b)-rib) 
(2.34) = &//(&). 

To find the critical points of fx, we differentiate: 

1 — 6 cot b -rj-r - cos 6- cos 6(1 - 6 cot 6) 

= 2 . 2 + a ^ rn-^ 

sm o sm o 

1 — 6 cot 6 6 sin 6—2 cos 6(1 — 6 cot 6) 

= 2 5 1" a o • 

sin 6 sin 6 

So 

1 — 6 cot 6 



(2.35) //(&) = ^ a 



cos 6(1 — 6 cot 6) — | sin 6 
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From this and (|2.33|1 we get: 

l(b m ) - b m n(b m ) = b m cotb m 

sm a 



= 1 - (1 - 6 m cot6 m ) 7-^7— a 

sin b m 

= 1 - a (cos b m (1 - b m cot b m ) - %^ sin b rn + ^ m ) 
V 2 sin b m J 

= 1 — a (cos b m (1 — b m cot 6 TO ) + sin b rn + 6 m cos b m cot 6„ 



V + ^ sin6 m ^ 



1 — a ycos 
1 - aS(b m ) 



For the proof of the next theorem we need another lemma: 



□ 



Lemma 2.6. On each interval (mir, (m + l)ir) , m G N 8 has a unique zero, and a 
unique critical point at b m , which is a maximum with S(b m ) > 0, if m is even and 
a minimum with S(b m ) < 0, if m is odd . 




Figure 4. The function 5. 



Proof. Since 



5(b) = <^ cos b + ^ sin b = 



<^ 2 + fetanb = 0, 
there is a unique zero in each interval (mir, (m + l)ir) , m £ Nq. 
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The critical points of S are given by 

5'(b) = - sin b + i sin b + ^ cos b = 
6 cot b - 1 = 
& b = 6 m . 

Since 

(2.36) 5" (6) = -^sin6, 

we see that & m is a maximum, if to is even and a minimum, if to is odd. □ 

Theorem 2.3. Given (xi,0), (x,u) € R n x R wi£/i xi ^ ±x and zt > 0, there are 
finitely many geodesies joining these two points. These geodesies are given by l|2.6p . 
where b is a solution of 

2u 

(2-37) a a =MW, 

Fi| + Pl 

and i/ieir lengths are strictly increasing with b. Moreover the shortest geodesic 
joining (xi,0) and (x, u) is given by the unique solution b € (0,7r) of (|2.37p in £/ie 
interval (0,7r). Wtt/i i/i«'s solution b the Carnot-Caratheodory distance is 



(2.38) dec ((zi,0), (a:, u)) = Lfr 6 ) 

_ b 
sin & 



xi\ 2 + \x\ 2 — 1x\ ■ xcosfr 



^2bu+ (|xi| 



2 



|x|" ) 6 cot 6 — 2 -xi • x. 

sin o 




Figure 5. Geodesies j 



loing different points. 
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Proof. From Lemma [2711 we have that the geodesies in this case are implicitly given 
by solutions of the equation -, — , y = [i(b). The properties of the function fi 

I *^ 1 I "T" | | 

show us that there is exactly one solution in the interval (0,7r) and at most two 
solutions in each interval (rmr, (m + 1)tt), to € N. So it remains to prove, that the 
values of I at these solutions are strictly increasing, i.e. if s < t are two solutions, 
then l(s) < l(t). 

Let s £ (to7T, (to + 1)tt) , to € N, then we have by (??) 



l(s) = I l'(b)db + l(b m ) 
bn'(b)db + l{b m ) 
bn(b)\i - f n(b)db + l(b m ) 



sfi(s)- fj,(b)db + l(b m ) - b m ^(b m ) 
s/j,(s) — / fi(b)db + 1 — aS(b m ). 



So, if we have two solutions s < t of (|2.37p in the same interval (rmr, (to + 1)tt) , to G 
N, we see that 

Z(s) - = sfi(s) - / /i(&)d6 - tfi(t) + [ n(b)db 



= {s-t)n(s) + / fi(b)db 

J S 

< o. 

Now let s m ,t m € (to7t, (to + 1)tt) , s r „ < t rn , to S N be the solutions of <|2.37[1 . To 
compare the corresponding "lengths" I at t m and s m +i, notice that 

(2.39) Z(t m ) - Z(»m+i) 

/■im rbm+1 

= - (s m +i - tm) (J.(t m ) - / fi(b)db- [i(b)db - (6{b m ) - 5(b m+ i)) a. 

Now consider the following cases: 

1. case a > and to even: 

Then, since mir < b m < b m , S(mir) = 1 and b m is a maximum of 6, it follows 
that 5(b m ) > 1. And, since & m +i < 6 m +i < (to + 2)7r, <5((to + 2)ir) = 1 and 
& TO +i is a minimum of 5, we have S(b m ) < 1. So (S(b m ) — 5(6 m +i)) a > 0, which 
gives the claim, since /j, > 0. 

2. case a < and to odd: 

Then, since rmr < b m < b m , 5(mir) = — 1 and b m is a minimum of 8, it follows 
that S(b m ) < —1. And, since b m +\ < b m+ \ < (to + 2)ir, S((m + 2)ir) = —1 and 
b m+ i is a maximum of S, we have 5{b m ) > — 1. So (5(b m ) — S(b m+ i))a > 0, 
which gives the claim, since fx > 0. 

3. case a > and to odd: 

Here 6 m < 26m--i < s m+ i < 6 OT +i. Let 

s . = /<m ,ifA(¥) >A(^) 

\ Sm+1 ,if < AC^) ' 
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Then, since 2b m -i is a minimum of jx (|) on (m7r, (m + 2)n) and monotone on 
the left and right of 2b m -i : max fc(E [ tm , Sm+1 ] fi (|) = /t (§) and therefore 

-(Sm+l - t m )n(t m ) = - (s m+ i - t m ) n(s) 

= - (sto+i - *m) ((1 - a)A(s) + a/i 

< - (s m +i - t m ) a/t (0 

r m+i - „ 

< —a / u - ) do. 
~ Jt m W 

Notice that has no singularity at (to + 1)tt. Now (|2.39[) becomes, since 

>a/i(f): 



Z(*m) - K s m+l) <—a I A ^2 J d& - (<*(&m) - 5(b m+1 )) a 

bm+1 ( 6/2 b\ 

I ^j+tan-J db - (6(b m ) - S(b m+1 )) a 

-a&tan^™ +1 - (S(b m ) - S(b m+ i))a 
a ( 5(b m+1 ) - 6 m+ itan^±i - ( S(b m ) - 6 m tan 6 



2 



Further 



5(b) — b tan - = cos b + — sin b — b tan - 

, • b ( b 1 A „ . 2 6 
= o sm - cos r — 2 sin - 

2 V 2 cos I / 2 

i ib b , 6\ 

- ( sm 2 - tan - + 2 sm 2 - J + 1 



2 sin' - 1 + - tan - 



-:„2 

2 

Since 6 TO < 2b m-i < 6 m+ i and 1 + L-i tanb m-i = and 6^1 + I tan If is 

2 2 2 Z Z 

increasing on (m7r, (to + 2)tt), we conclude: 

l(tm) ~ l( s m+l) 

. n ( ■ 2 °m ( 1 , b m b m \ . , 6 m +i / 6 m +l &m+l 

< 2a sm — 1 H tan — — sm 1 H tan 

V 2 V 2 2 J 2 \ 2 2 

< 0. 

4. case a < and to even: 

This is very similar to the third case, if we use ft. Here b m < 2b™. < s m+ i < 
b m+ i. Let 



s := 



t m , if fl^-f) >£(*^) 
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Then, since 26™ is a minimum of fi (|) on (rmr, (m + 2)n) and monotone on 
the left and right of 26™.: niax;, e [t m)STO+1 ] = /5(f) and therefore 

-(%+i - t m )n(t m ) = - (s m+ i ~ t m ) fj,(s) 

= - (sm+i - t m ) ((1 + a)fl(s) - afl 

< {s m +l - t m ) Cljl ^-^ 

f-Sm+1 / b \ 

< a Jl I — ] do. 



Notice that has no singularity at (m + l)ir. Now <|2.39|) becomes, since 

fi(b) > -aft(b/2): 

l{t m ) - l(s m+ i) < a ^ fi db - 5(b m ) + c5(6 m+ i)^ 

hm+1 ( b/2 b\ \ 

-r^Tb - cot - tft - 6(b m ) + S(b m+1 ) 



a ( -6cot^|£™ +1 - 8(b m ) + 8(b m+1 ) 



Further 



" ( b m cot ^ - 8(b m ) ~ (b m+ i cot - 5{b m+ i)) 



b cot 5(b) = — cos 6 sin 6 + 6 cot - 

2 w 2 2 



6/6 1 \ 2 6 

-6 cos - | sin 7- I — 2 cos — + 1 

sin - 



2 I 2 b 



-2 cos- - f L - -col - I -! I. 

b „„ + b 



Since 6 m < 26221 < 6 m+ i and 1 — 6222 cot 62= = and & ^ 1- - cot f is increasing 
on (m7r, (m + 2)71), we conclude: 

l(trn) — l{ s m+l) 

. n / 2°m ( -, b m b m \ 2 6 m +i / 6 m+ i 6 m +i 
< —2a cos — 1 cot — — cos 1 cot 



2 



0. 



If <o € (0, 7r) denotes the solution of l|2.37p . it remains to show that l(t ) < l(si). 
The above calculations can be extended to this, if one sets 60 = 0. In detail we 
have 

/(to) = f° l'(b)db + l(0) 
Jo 

/■to 

= toM(*o) - / n(b)db + 1 - a, 
Jo 

so that 

l(to) - l(si) = -(si - io)M(si) - / ° Kb)db - I ' ^{b)db - (1 - 5(6!)) a. 
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Again from the previous lemma one gets l(to) < l(s±) if a > 0, since l — 5(bi) > 
and /j, > 0. If a < 0, then since fi is increasing on [0, 7r): 

-{si - to)j"(si) = -(si-to) ((l + a)A(si)-aA(yJJ 
(si - i )a/i (yj 



< 



So we get 



< a j fi[-)db. 



2 
b 



l(*o)-I(*i) < a/ A ( - Jtft - (1 - 5Qn)) a 



a ( -focot^ 1 - 1 + S(bi) 



a [ Mh t j 6iCOty + 1 

= 2a cos — 1 — — cot — 
2 2 



< 0, 

since 1 - § cot \ > on [0, 2vr). 



□ 



Theorem 2.4. Given (x, 0), (x, u) € R n xM, i^0, u > 0, i/iere are finitely many 
geodesies connecting [x,0), (x,u). Namely: 

1. For any solution b of 



(2.40, 

there is a geodesic 7 b given by (|2.6p with square of length 

(2.41) L{ n b f = \x\ 2 = 2bu~2 \x\ 2 fetan 

cos" 1 1 2 

which strictly increases with b. 
2. If u > \x\ 2 tt, then for each m E N with u > mir\x\ 2 there are geodesies 7 b,c 



given by (|2.6p wii/i 6 = 2mir, \c\ = 2\J mnu — (ran) 2 \x\ 2 , with lengths 
(2.42) Hl b,c ) = 2^fmnu. 

Moreover the shortest geodesic and therefore the Carnot-Caratheodory distance is 
given by the unique solution b of (|2.40p in the interval (0,7r), with this b we have 



(2.43) d CC ({x,0),(x,u)) = ~^-j\x\ = \ 2bu - 2 \x\ 2 fctan 



cos ■ 



Proof. Notice that in this case a = 1. We will consider two cases: 
1. case b 6 7rN: 

Then since 7m (1) = x > there exists to £ N with 6 = 2m7r. Then 7,^ (1) = u is 
equivalent to 



(|c| 2 + (2m^N 2 ) i i 7r . 
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1 2 3 



Figure 6. Geodesies joining (1, 0) and (1, 6) with different b g Z* 
vs. the two geodesies with b = 2n and c = ±y/2bu — x\b 2 . 

So if u > mir \x\ 2 , there is a geodesic given by (|2.6p with 
6 = 2to7t, |c| = 2\/ 7717™ — (to7t) 2 |x| 2 , 

with length: 

L 2 ( 7 b < c ) = |c| 2 + \x\ 2 b 2 = Amnu 

This gives the second part of the theorem. 

case b $ ttN: 

In this case we have 

Kb) = A (J) ■ 

This is very similar to the situation in Theorem 12. 3} but much easier, since we 
are only dealing with /i. From Lemma l2~31 we know that the equation 




has at most two solutions in an interval ((2m + l)ir, (2m + 3)?t) , m £ No, each 
defining a geodesic. If s is a solution of this equation, the square of the length 
of the corresponding geodesic 7 is given by 

l 2 ( 7 ) = 2;( s )N 2 

s 2 2 
= 2 ' „ (1 — coss) |a;| 
sin s 
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If s e ((2m + 1)tt, (2m + 3)n), m e N , we have 
l(s) = { l'(b)db + I (2bjj 



2b„„ 

bfi' (b)db + I (2b r , 

>-b m \ 

b -fi' (J ) db + l [2b„ 

2b„ 



2^ \2 



(2.44) = api (|) - 2b m pi (& m ) - fit db + I (2b n 
But now 

2b 2 - ( b ~ \ 

l(2b m ) - 2b m (i(b m ) = - 2b m ^— + tan& m 

cos^ b m y cos^ o m / 

2 ^77^ t an bffi 

(2.45) = 2. 

First let s < t be two solutions in an interval ((2m + l)n, (2m + 3)tt) , m e N . 



Then by (??) and (??) 

- I(t) = «/!(-) - I ,,(-),//,-/,,(-) + / /, ( - j ,lb 



2b m \^/ V Z / J 26 



< o. 

Now let t m € ((2m + 1)tt, (2m + 3)n) and s m+ i € ((2m + 3)tt, (2m + 5)?r) be 
two solutions, where t m is the rightmost and s m+ i is the leftmost solution in 
the corresponding interval. Then by (??) and (??) 

l(t m ) — l(s m+ i) 



2b m V z / \ Z / J 2b . \Z 



< o, 

since 26 m < i m and s m+ i < 2b m+i . 

If t € (0,7r) is the solution of (|2.40p . we have to check, that l(t) < 1(sq). As 
above, we have 

l(t) = [ l'(b)db + l(0) 

= tf*(z) - I v(^)db + 2, 



so that 

no -w (|) -jf a (|) *- «»A (f) + (| 



/ \ f\ fb\ „ /" 2bo „ / /« 

so 



< 0. 
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At last we have to compare the length of the shortest geodesic of case 1 (if there 
is one) with that of case 2. Let b G (0,n) be the solution of (|2.40p . then since 
tan | > 0: 

L 2 ( 7 b ) = 26w-2|a;| 2 6tan^ 

< 2wu 

< L 2 ( 7 2,r < c ). 

□ 

Theorem 2.5. Given (x, 0), (— x, u) G M™ x M, x ^ 0, u > 0, there are finitely 
many geodesies connecting (x,0), (— x,u). Namely: 

1. For any solution b of 

< 2M> Wit 

there is a geodesic 7 b given by (|2,6p with square of length 



(2.47) L( 7 b ) 2 = \x\ 2 = 2bu + \x\ 2 6cot ~ 

sm § 2 



which strictly increases with b. 



2. If2u > \x\ 2 tt, then for each m € No with 2u > (2m+l)7r \x\ z there are geodesies 



7 b < c given by (HU with b = (2m + l)7r, |c| = \j2(2m + 1)ttu - (2m + 1) 2 tt 2 \x\ 2 
and lengths 



(2.48) L{^ c ) = v/2(2m + 1)ttm. 

Moreover the shortest geodesic and therefore the Carnot-Caratheodory distance is 
given by the unique solution b of (|2.46p in the interval (0,7r), if 2u < \x\ 2 ir. With 
this b we have: 



b , L, , 2 , .6 



(2.49) d CC {(x,0),(-x,u)) = -r-^N = \/26u+|a;r&cot- 



sm 



2 



v4nd, i/2u > |x| 2 7r, then the Carnot-Caratheodory distance is given by 

dec ({x, 0), (-x, u)) = \J2-ku. 



Proof. Notice that in this case a = — 1. The first part is the same proof as above, 
if one replaces fi, b with /i, b. In particular: 

1. case 6 G 7rN: 

Then since 7(i)(l) = x, there exists m G No with b — (2m + Then 
7(2) (1) = u is equivalent to 

So if 2u > (2m+l)7r |a;| 2 , there are geodesies 7 b,c given by l|2.6p with b = (2m + 
|c| = ^2(2m + l)nu - (2m + 1) 2 tt 2 \x\ 2 , with length 

L 2 (^ c ) = \c\ 2 + \x\ 2 b 2 = 2(2m + l)%u 
This gives the second part of the theorem. 
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2. case b $ 7rN: 

In this case we have: 

fi(b) = ~J\ 

This is again very similar to the situation in Theorem I2.3[ but a lot simpler, 
since we are only dealing with jl. From Lemma [2731 we know that the equation 



(2.50) = ft 



has at most two solutions in an interval (2to7t, (2to + 2)7t), m £ N, each defining 
an geodesic. If s is a solution of this equation, the square of the length of the 
corresponding geodesic 7 is given by 



l 2 ( 7 ) = 2;( s ) \xY 



sin 2 s 



(1 + coss) \x\ 



sin 2 f 



\x\ 2 . 



If s e (2m7r, (2m + 2)ir), m G N, we have: 
l(s) = I l'(b)db + I (26 m ) 



2b„„ 

bfj,'(b)db + 1 ( 2b 

2b m ^ 

b_, (b 



l2b m 2 
Sfl 



y! f^jdb + l(2b m ) 
(0 - 2b m fi (& m ) - ^_ db + I (2b„ 
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But now 

l{2b m ) - 2b m jl(b m ) = -^T- ~ I — TT ~ cot ~ h ™ I 
sin b m \sm b m J 

= 2b m cot b m 
= 2. 

First let s < £ be two solutions of (|2.50|) in an interval (2m7r, (2m+2)7r), m g N; 
then 

i W - iw - •/* (0 - £ a Q) * - m (0 + £ a Q) * 
= -(*-*)a(~H f a(Tw& 



< 0. 



Now let t m S (2to7t, (2m + 2)7r) and s m +i € ((2m + 2)tt, (2m + 4)tt) be two 
solutions, then 



Z(*m) - l( s rn+l) 
&m A 



t) - jC A G) d6 - sm+iA + i-J! a G) * 



2 / -/26 m \2/ i Sra + 1 V - 



o, 



since 26 m < i m and s m +i < 26 m . 

If to € (0, 27r) is the solution in the interval (0,27r), we have to show that 
l(to) < l(s%). But the above calculations show that 

l(t )= f°l'(b)db + l(0) 



to 



to 



*oA ( ~ ) ~ I A ( 2 ) d6 + 2 - 



So we get 



'Co) /--m ) = /„// ( |) - ^° A Q) A ~ «iA (y) + £ A (T) 



2 J ./ n f V 2 y ,/ Sl " V 2 



< 0. 



At last we have to compare the length of the shortest geodesic 7 7r '°, |c| = y 27ra — 7r 2 |x|' 
of case 1 (if there is one) with that of case 2. Let b £ (0, 2ir) be the solution of 
(|2.46p . If b < 7T, then, since A is monotone increasing on (0,7r), we get: 

FF = "U 
< >' (I 

7T 
= 2' 
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But this means, that 2u < ir \x\ 2 and therefore there is no extra geodesic. If now 
b > it, we have to show that L 2 (-f b ) > 2iru. We have 



u _ t b 
2 

b 



2 



b 

- cot- 



2 & 



cot - [ ^ 1 

CO 2 I 2 sin | cos t 



2 2 



so that 



6 — sin b b 

— ^~T~ cot o' 

smo 2 



L 2 ( 7 6 ) = 2ou + 2|a;| 2 ocot^ 
sin 6 

2bu 1 + 



2u- 



6 — sin 6 
& 2 



6 — sin 6 
> 2tto. 

The last inequality holds since /(&) := & _g in b > 7r, for all 6 <E [71- , 2n]. To see this, 
observe, that f(ir) = ir and / is monotone increasing: 

, _ 2b{b - sin b) -b 2 (l- cos b) 
■' { } ~ (o-sino) 2 

_ o 2 (l + cos6) - 2&sin6 
(b — sino) 2 

> 0. 

This completes the proof, since the geodesies j b ' c in case 2 (with parameter b = it, 
\c\ = \J2ku — 7r 2 x 2 ) are shorter. □ 

To close this section, we give a summarizing theorem about our results of this 
section: 

Theorem 2.6. Let (x\,u\), (x,u) G R™ +1 . If x\ = —x and2\u~ u\\ > tt \x\ 2 , 
then the Carnot-Caratheodory distance is given by 



(2.51) dec ((xi,U!), (x, u)) = ^2n \u - m\. 

Otherwise the Carnot-Caratheodory distance is given by 



b I 2 2 

(2.52a) dec ((x\,u\),(x,u)) = - — rvl^il + M — 2xi-xcos6 

smo v 



(2.52b) = \/2b(u-u 1 ) + (|xi| 2 + |x| 2 )6 cot b - 2x 1 ■ x-^— 

V smo 

where 6 e (— n, w) is the unique solution of 

2(u — ui) b x\ ■ x l — bcotb 

2.53 -\ 2 = -^--cotb + 2 1 — — 

\xi\ +\x\ sm 2 6 \ Xl \ +\x\ sm& 

in the interval (— 7r, it). 
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3. Heat kernel 

We now turn our attention to the heat equation, i.e. the Cauchy problem 



(3.1) 



(d t - G)v(t,x,u) = 
v(0,x,u) = f(x,u) 



with some function / £ iS(R™ +1 ). Taking the Fourier-transform in the it-variable 
we see that this is equivalent to 



(3.2) 



(d t - H\)w(t, x, A) = 
w(0,x,X) = g(x,X) 



where g{x, X) := f{x, X) denotes the partial Fourier-transform of / in the u-variable 
and H\ — A x — |A| 2 |a;| 2 denote the Hermite-operator with parameter A £ R, which is 
the partial Fourier-transform of G. Thanks to the Plancherel formula, any solution 
w of p.2p defines a solution v to l|3.ip by taking the inverse Fourier-transform, and 
vice versa. 

But the latter equation is well known and we give a brief discussion of this 
equation. The solution of this equation is given by e tHx g, where e tHx is a bounded 
operator defined by the functional calculus, since H\ is a self-adjoint operator. By 
the Schwartz kernel theorem e tHx is given by an integral kernel, i.e. 

(3.3) e tH »g(x)= [ fc t A (x, Off(0«- 

And, since H\ is hypoelliptic, we have k$ £ C°°(R 2n ). 
The Hermite-functions 

nn/4 / _l \\ a \ / r) \ a 

(3.4) h a (x) := e-M (-) e^-l , x £ R", a £ NJf, 

which give an orthonormal basis of L 2 (R"), are eigenfunctions of H 2 -„: 

(3.5) H 27r h a = -2n(2\a\ + n) h a . 
Since 

(3.6) (A - M 2 |x| 2 ) f(r-) = - 2 ((A - (^) 2 | • | 2 ) /) (r.), 
we can rescale the Hermite functions 

to get an orthonormal basis (ha) ae ^ of eigenfunctions of H\: 
(3.8) H x hl = -\X\{2\a\+n)h x a . 
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Using Mehler's formula, we can calculate the kernel & A : 

\a\>0 



|a|>0 



tt(1 + e" 4 ^ 4 ) ^™|A| + 4 7 rf!|A| e - 2 l A l t> 



• exp 

(27r)-"/ 2 



1 - e — 4| A|* 

|A| 



sinh(2|A|t) 



■exp(-^^|A|coth(2|A W + -^,.e N 



(Ant) 



-n/2 



2\t y 

sinh(2Ai) J 



•^.(~(2A*oaa(2At)(|x|^|e| a )- s ^*.f 



Now, by taking the Fourier inverse 



(3.9K f (i,^) = (2i)- 1 / fc t A (a;,Oe lA "dA 



= (27r) -1 (47rt) 



/ 2\t \ f 
^sinh(2At) J 



^ <-x P ( -~ (|2Aicoth(2At)(| a; | 2 + |e| 2 ) - — ) ) < lX " <IX 



i/2-l f ( A 



\ sinh A 



exp | ^AcothA(|a;| 2 + |£| 2 ) - ■ £- 2iAu j ) ./A. 



one gets the heat kernel for the Grusin operator. Now the solution of l|3.ip is given 

by 



(3.10) v(x,u)= / K t (x,Z,u-\)f(Z,\)dZd\ 
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4. Estimates 

We may now use the results of section [2] to give some estimates for the heat 
kernel. We set 
(4.1) 

Kx,t,u) := J f-A^Vexpf- (^AcothA(|x| 2 + |C| 2 )--^^-e-2*A U ^dA 

— OO 

OO n 

^ A ' exp(-(AcothA — — -a J R 2 + 2iXu ) dX 



with 



\ sinh X J \ \ sinh A 



/? : ^ /?(.r. i ) := V \x\ 2 + |£| 2 and a := a(x, £) := ^ G [0, 1] . 



A 



Then 

(4.2) 

It is also convenient to set 

^ ^ ' 'sinhA 

(4.4) V(A) := ^(A,o) 

:= AcothA r-T-r a ' A G K. 

sinh A 

One may see that ip{ib) — bcotb — j-^bja, so that if b is a solution of (|2.53p . then 
the exponent 

A coth A t— ra ) R 2 - 2iAw 

sinh A / 

at A = ib gives exactly the square of the Carnot-Caratheodory distance dcc((x, 0), (£, m)) 2 . 
So we expect Gaussian-type estimates of the form 

where F(x,£,u) > is a function depending on x, £,u. It will turn out that F has 
polynomial growth, more precisely 

F(x,^u)<(l + dcc((x 1 0),(^u)fr, 

where a = % — 1, if n > 2 and a — if n < 2. Compared to the euclidean case, 
where a = 0, one has some additional growth, if the dimension is greater than 2. 

The strategy is as follows: Move the line of integration in l|4.ip from the real 
axis to the line R + ib, where b is near to the solution bo of (|2.53p . 

Lemma 4.1. The function 

v ^ = {Ax 

where the square root is the principal branch in C\ (— oo,0], is holomorphic in 
{z G C : \^sz\ < ir}. And for v + ib G C, — n < b < n, one has 

(4.5) \V(v + ib)\ < 



sin b , 

b 2\~> 
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Proof. To show that V is holomorphic, we have to show that 

for all A e C with -tt < 3?A < tt. Let z/ + ib e C, |6| < tt. Then 
v + ib (v + ib) sinh(^ — ib) 



(4.7) 
So 



sinh(V + ib) | sinh(i^ + z6) | 2 

v sinh v cos b + b cosh v sin 6 + i(b sinh z> cos b — v cosh i/ sin b) 



i/ + ib 



| sinh(f + ib)\' 



€ R <^4> 6 sinh z/ cos 6 = v cosh zv sin 6 



sinh(z/ + i&) 

6 cot b = v coth z/ V v = 0. 

But b cot 6 < 1 and z/cothz/ > 1 for v ^ 0; and f = means 

v + ib 



sinh(z^ + ib) sin b 



> 0. 



z/ + ib 



We can conclude 

( - 4 ' 8 ' 1 sinh(z^ + ib) 

The second inequality is easy, since 

v + ib 



— ^(-oo,0], v e R, -tt < b < TT. 



sinh(z/ + ib) 



sinh 2 v + sin 2 b 



< 



sinh 2 v 
b 2 



= 1 + ^7 



/ V 



v 2 ] \ sinh v 



To verify the first inequality we show that 

...... 2 

m ■■= 



v + ib 



< 



sin 2 b 



, for all v <G 



sinh(z/ + ib) 

This is obviously true for v = 0, and since f(—v) = ,f(v), we may restrict our 
analysis to the case v > 0: 



/» 



2z^(sinh 2 i/ + sin 2 b) - (y 2 + b 2 ) sinh(2z^) 

(sinh 2 v + sin 2 b) 2 
2v sinh z/(sinh v — v cosh z/) + 2v sin 2 b—b 2 sinh(2z^) 



(sinh^ v + sin 2 6) 2 



< 0, 



since sinhf < z^coshz^ and 2z^sin 2 6 < sinh(2z^) sin 2 b < 6 2 sinh(2z/). This means 
that / stays below for all v e R. 

Lemma 4.2. T/ze junction 



□ 



(4.9a) 
(4.9b) 
(4.9c) 



tp(X) = A coth A — a 



A 



sinh A 



(1 — a)AcothA + aAtanh 
(1 + a)AcothA — aAcoth 



HEAT KERNEL ESTIMATES FOR THE GRUSIN OPERATOR 



28 



is holomorphic in C\ {irki : k 6 Z*} for a 7^ ±1. 

If a = 1, £/zen i/>(A) = A tanh -| is holomorphic in C \ {(2fc + l)-7ri : fc G Z}; 
if a= — 1, i/ien ?/>(A) = A coth -| is holomorphic inC \ {2k-xi : fc € Z*}. 

In particular if) is holomorphic in the strip {z 6 C : < 7r} /or an?/ — 1 < a < 1. 

Proof. The holomorphic properties of -0 follow easily from those of the trigonometric 
function coth, tanh, and the equalities l|4.9p . which are easy to verify: 

tp(X) = A coth A — a— ^ 



and 



sinh A 

/, x x , x > cosh A — 1 

= (1 — a) A coth A + aA — — - — 

sinh A 

= (1 — a)A coth A + aA tanh — , 



, . „ , . cosh A + 1 

0(A) = (1 + a) A coth A - aA — — — 

sinh A 

= (1 + a)A coth A — aA coth — . 



Lemma 4.3. For v + ib € C, |6| < tt, ip(ib) is real and 
(4.10) $l(ip(v + ib)) > ip(ib), 

where 

(4.11a) ip(ib) = bcotb — a 



sin b 

b 



(4.11b) = (l-a)6coto + a (^-6tan 2 

(4.11c) = (1 + a)6cot b — aocot ^. 

Proof. dHU]) is clear. Let v + ib e C, |b| < tt, then 

«.,, .,s , ., NN 1 z^smh(2z/) +6sin(2o) 

3? {{v + ib) coth(z/ + ib)) = , 2 ; 

2 sinh v + sin b 

v coth v sinh 2 z/ + & cot b sin 2 6 

sinh 2 v + sin 2 6 

2 v coth v — b cot 6 

= cot + sinh 2 5 — 

sinh + sin b 

> b cot b, 

and for |6| < f : 

1 i/sinh(2i>) - 6sin(2&) 



S? ((z/ + i&) tanh(z/ + = 



2 sinh 2 z/ + cos 2 b 
v coth z> sinh 2 v — 6 tan 6 cos 2 
sinh 2 z/ + cos 2 b 

. 2 v coth z^ + 6 tan 6 
— otano + smh v- 



sinh 2 z> + cos 2 b 
> —6 tan b. 



□ 
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This gives for A = v + ib £ C, \b\ < tt: 
K(V(A)) 



3?((l-a)AcothA + aAtanh|) , if a > 
3? ((1 + a) A coth A - aX coth |) , if a < 



> 



(l-o)6cot6 + o(-6tan|) , if a > 



, (1 + a)6cot6 — afecot | 
^»(i6). 



, if a < 



□ 



With these preliminary lemmata we can start estimating h by moving the path 
of integration to A i— » A + ib: 

Lemma 4.4. There is a constant C > 0, such that for all x,£ £ R™, u £ R and 

6 £ WL, \b\ < tt the following estimate holds: 

(4.12) 



\h(x,£,u)\<C 



exp — 2bu — bcotb- 



sin b 



-co 

CO 



sin b 

with a = max (0, J — l) . 
Proof. We have 

GO 

y \V{\ + ib)\e-^ x+ ^ R2 - 2bu d\ 

— CO 

CO 

= e -^( 4 6)fl 2 -2b„ J \ V (X + ib)\dX. 

— CO 

Now we split the remaining integral with r := it — \b\: 

CO CO 

J \V(X + ib)\dX = 2 J \V(X + ib)\d\ 

-oo 

r n 7r Ti oo 

*/i^r-/('4)'-/ 



(kl 2 + KI 2 ) , 



< 



< 



< 



< 



sin b 
b 

sin 6 
6 



A 



sinh A 



dX 



A™/ 2 



rfA + 1 



sin& 
sin& 



+ 1 
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since 

6 1 1 



(4.13) 



sin 6 tt — \b\ r 

□ 



Observe that we can assume u > 0, since h(x, £, —u) = h(x, £, u). 

Setting b — bo, where bo parametrizes the corresponding (shortest) geodesic, 
would not give good results in every case. If b tends to 7T, the estimate would blow 
up. To avoid this problem, just set b = bo — e with an appropriate e > 0. (Since we 
assumed u > 0, we have b > 0. 

Proposition 4.1. For £ = (x, 0), r] — (£,lt) G R n+1 we have 

(4.14) \h(x, u)| < min (l + ^ C | C ^ ?) , 1 + d cc (C, r?) 2 ) " e~ d °°^\ 

with a = max(^ — 1,0). 

Proof. First assume that x ^ — £. Then the shortest geodesic and therefore the 
Carnot-Caratheodory distance will be parametrized by o G (7r,7r). We claim that 

h l <3 ( 1 , dcc(C,V?\ 



sin 2 6 y \x + £\ 2 J ' 
To see this, first assume that a < 0, i.e. £ • x < 0. Then 

dcc(C, V) 2 = -4- (N 2 + ICI 2 - 2^ • x cos bo) 
sm 6 v ' 

>-^-(N 2 + iei 2 + 2C- 

sm Oo v 
h 2 

'^ + e| 2 . 



sin 2 6o 

If a > 0, i.e. £ • a; > 0, then 



dcc(Cv) 2 > \-%-(l-acosbo)(\x\ 2 + 2t-x + \t\ 2 
o sm Oq v 



1 / 6q OqCOSOq 



o \sm Oo sm oq 



2 



If 6 > f , then cos6 < and if b > § , then cos6 < < \ 

So 



1 f bl 



dcc^) 2 > L A^--'-\V + t? 
o \ sin On 



1 f bl 



.?--3 i.+ei 2 . 

o \sm Oo / 

This proves the claim. Now use Lemma [HU with 6 = 6 to get 



(4.15) 



\h(x^u)\<^ + ^f^ye-^y 



But this remains true, if x = — £ (and x = £ = 0), since the right side becomes +oo 
in this case. To prove the remaining estimate 

(4.16) \h(x,£,u)\ <(l + d C c(C^) 2 ) a e- dcc ^ 2 , 

we, again, use lemma 14.41 Let bo G [0, n] be the parameter for the Carnot- 
Caratheodory distance. As mentioned above assume that u > 0. If bo < § , then 
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si ^° feo < -|, and one can immediately use Lemma \4. 41 with b — bo. So assume that 
bo > § • Now set b = bo — e with e < f , which will be specified later. 
If 6q ^ we have 

b cot b = bo cot 60 + 6 (cot b — cot 60) — e cot 60 

> 60 c °t bo — e cot 60, 

since cot is monotone decreasing on (0, 7r); and 

6 bo , f bo b\ bo 

—0 tan — = — On tan ho tan tan — + e tan — 

2 2V 2 2 J 2 

bo bo 

> —On tan he tan — , 

- u 2 2 ' 

since tan is monotone increasing on (Q,7r). Therefore 



(1 — a)6cot6 + a(— 6tan|) , a>0 
(1 + a)6 cot b — ab cot | , a < 



oo 



The use of lemma I4~41 with this b yields: 

(4.17) MX,£,U)\ < (-^) Q e -Wbo)R 2 -2b u +t U,(ib ) + 2 bo u) 



first for bo 7^ 7T. But, if 60 = tt, then we have x = £ = 0orx = — £. We claim that 
in this case estimate l|4.17p remains true: If x = £ = then 



\h(x^, U )\ < ( _±_) a e -2boU+2b u^ 
\sm b J 



_ 1 e -(l-^)dcc(f ;t 7) 2 



since 26nu = 27ru = dcc{C, V) 2 ■ And, if x = — £, then a = — 1 and i/i(io) = ocot | > 
0, so that 

1, / ^ \\ ^ ( t> \ -4>(ib)R 2 -2b ii+2b u^- 

\h(x,£,u)\ < — - e 
\ sin / 

< e -26 M+26o« T 

1 

c 

since again 26nii = 27ru = dcc{C, V) 2 ■ 

Setting e = ppg~7g t? )a , which optimizes this inequality (up to a constant) , gives 
us the desired result: 

(4.18) \h(x,£,u)\ < (l + d C c(C,v) 2 re- d ^^ 2 . 

□ 

This gives the following main result: 
Theorem 4.1. For £ = (x, 0), 77 = (£, u) € M™ +1 we have 

(4.19) € , U )i < t-t-i min ^ + *eM , 1 + ^|^) Q e -^-(^ , 
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with a = max — 1,0). 

Proof. The proof is almost trivial, if one uses (|4.2j) and the homogeneity of dec ■ O 
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